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Abstract — This paper focus on design of multivariable 
controller for Quadruple Tank Process, a two input two 
output system with large plant uncertainty using QFT 
methodology. In the present work, a new approach using 
Quantitative Feedback Theory (QFT) is formulated for 
design of a robust two degree of freedom controller for 
Quadruple Tank Process. The design is done in frequency 
domain. This paper presents a design method for a 2 x 2 
multiple input multiple output system. The plant 
uncertainties are transformed into equivalent external 
disturbance sets, and the design problem becomes one of 
the external disturbance attenuation. The objective is to 
find compensator functions which guarantee that the 
system performance bounds are satisfied over the range 
of plant uncertainty. The methodology is successfully 
applied to design a two degree of freedom compensator 
Quadruple Tank Process. 

Kewords — Quadruple Tank Process, Robust Control, 
QFT, MIMO, EDA. 



I. 



INTRODUCTION 



Many practical applications are multi- input multi- 
output (MIMO) systems. MIMO system is composed 
of several interacting control loops and numbers of 
alternative configurations of control loops are very 
large. Interaction produces several undesirable effects. 
The design method proposed in this paper consists of 
design of multivariable controller to have desired 
response from the plant in the presence of plant 
uncertainties and external disturbances. The quadruple 
tank process is a multivariable laboratory process 
developed by K.H. Johansson[7]. This process consists 
of four interconnected water tanks and two pumps as 
shown in figure 1. Here objective to control the level in 
the lower two tank with two pump. 

The Quantitative Feedback Theory (QFT) is a 
control method in frequency domain which uses the 
Nichols chart [1]. It is an engineering science devoted 
to the design problem with quantitative bounds on the 
plant parameters and quantitative tolerances on the 
acceptable closed-loop system response. The objective 
is to find compensator functions which guarantee that 
the system performance bounds are satisfied over the 
range of plant uncertainty. 



*2 




Figure 1 : Quadruple tank system. 

In this work, we propose a QFT based methodology 
to design a robust two degree of freedom MIMO 
compensator for quadruple tank process in the 
presence of parametric uncertainty. The plant 
uncertainties are transformed into equivalent external 
disturbance sets and thus the design problem becomes 
one of external disturbance attenuation (EDA). The 
method has been proven rigorous by means of 
Schauder's fixed point theorem [2]. The design 
technique is successfully applied to design a two 
degree of freedom compensator for quadruple tank 
process. 

The present paper is organized as follows: Section II 
deals with the problem definition for a 2 x 2 MIMO 
design problem. Basic concepts are discussed in 
section III. Design equations used in the synthesis of 
controller for quadruple tank process is given in 
section IV. Section V deals with the design procedure 
to design controller and prefilter. In section VI, the 
methodology is applied to design a compensator for 
quadruple tank process. Conclusions are drawn in 
section VII. 

II. A 2X2 MIMO DESIGN PROBLEM 

The two degree of freedom controller structure used 
for the design of two input and two output quadruple 
tank process is shown in figure 2. 
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Figure 2: Controller Structure for 2 x 2 MIMO System. 

P= [pij (s)] represents the set of all possible plant 
transfer function elements. If t(s) represents the closed 
loop transfer function matrix between Y(s) and U(s) 
and T(s) represents the overall closed loop transfer 
function matrix between Y(s) and R(s), then, 
T(s) = t(s)xF(s) 

Where, 
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F(s) is the forward filter designed to achieve the 
desired tracking properties for T(s). The controller 
placed inside the feedback loop has a diagonal 
structure given by, 



G = 



Six 

o 



Both G(s) and F(s) are rational and strictly proper 
linear time invariant matrices. 

III. BASIC CONCEPTS 

Consider a two degree- of- freedom controller 
structure as illustrated in figure 3, where P represents 
the set of transfer function which describe the region of 
plant uncertainty, G is the compensator and F is an 
input pre-filter transfer function. The system has two 
inputs: R the desired input signal to be tracked and D 
an external disturbance input signal which is to be 
attenuated to have minimal effect on output Y. 

The controller and prefilter are to be synthesized in 
order to meet robust stability and closed loop 
specifications. Compensator is designed so that the 
variations of output due to the uncertainty in the plant 
are within allowable tolerance and effects of 
disturbances are small. The compensator G is designed 
via loop shaping process. The 
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Figure 3: Two degree of freedom controller. 

loop shaping is done to generate a controller G(s), 
such that the open loop transfer function L(s) = GP 
satisfies certain specifications such as Nyquist stability, 
gain-margin, phase-margin etc. Once the compensator 
is designed the closed loop Bode plots will be as close 



together as the tracking specifications but they will not 
usually fall within the tracking specifications. A 
prefilter F is designed which brings the frequency 
response of the closed loop system within the tracking 
specifications. In figure 4, P is the set (due to 
uncertainty) of the n x n plant transfer function matrix 
P = [pij(s)]. The elements of the n x n compensator G 
= [gij(s)] and the n x n pre-filter F = [fij(s)] are to be 
chosen practically (each with an excess of poles over 
zero). They must ensure that the closed loop response 
falls inside the desired bounds over all P C P. An 
alternative method presented in [2] involved the 
replacement of the plant uncertainty entirely by an 
equivalent disturbance set and normal fixed, certain 
plant. 



F ^Q > G > p 

-I 



Figure 4: A Canonical structure. 

Consider a two- degree of freedom, n x n input- 
output feedback system, as shown in figure 4. Let P , P 
be the n x n nominal and perturbed plant transfer 
matrix, respectively. L = PoG is the nominal loop 
transmission matrix T = [t y o], T = [ty] are the nominal 
and perturbed system transfer matrix, respectively. The 
n x 1 output matrix c can be expressed as, 

c =PG{Fr —c) r^\ 

Where r an n x 1 input matrix, F is an n x n fixed 
pre-filter matrix of transfer function. 

(/ -hPG)c =PGFr (A\ 

P must be nonsingular matrix for the system to be 
controllable, so 

{p^ +g)c =gf 
{p p^ +p g)t =p gf 
{p p^ +l )t =l f 



Define 
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The nominal system transfer matrix 
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Figure 5: Disturbance attenuation only. 
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Equation 10 is implement in figure 5 with the 
equivalent disturbance matrix D e = VT. The uncertainty 
in P is transformed into the equivalent disturbance term 
VT; hence, it is called the "equivalent disturbance 
method". As T roams over its acceptable set, D e 
describes a corresponding set. 

IV. DESIGN EQUATIONS 

From equation 10, 
If L is taken as diagonal then 



At, 



V ij t i j + Y, v ik t k j 
k*i 
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-, 7 = U,3. 
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From equations 1 1 the design equations are derived 
as follows, For i = j 
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Substitute 
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Similarly from equation 1 1, for i ^ j 
1 
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where, 
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The normal loop transmission function ln should be 
designed to satisfy 
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(16) 



The Nichols chart with its loci of constant \L/(l + L)\ 
etc an be for finding the bounds on l iio to guarantee the 
desired performance. The bounds of lu (s) are 
calculated by using equation 15 and 16. Loop shaping 
is done such as at each frequency luo must lie on or 
above the specified frequency bound. 

V. CONTROLLER SYNTHESIS PROCEDURE 

The design procedure to design the prefilter F(s) and 
the compensator G(s), to achieve specific robust design 
for given region of plant parameter uncertainty is as 
follows : 

1 . Translate of time domain specifications into 
frequency domain specifications. 

2. Stability bound generation. 

3. Conversion of tracking specifications into 
disturbance. 

4. Disturbance bound generation. 

5. Grouping of bounds. 

6. Intersection of bounds. 

7. Design of controller G(s). 

8. Design prefilter F(s) for tracking of reference. 



VI. QUADRUPLE TANK APPLICATION 

The quadruple tank system is as shown in figure 1. 
Here objective to control the level in the lower two 
tank with two pump. The process inputs are vi and V2 
(input voltages to the pump) and the outputs are yi and 
V2 (voltages from level measurement device). The 
positions of the valves determine the location of 
multivariable zero for the linearized model. The zero 
can be put in either the left or the right plane. Thus 
system may have minimum phase characteristics or 
non minimum phase characteristics. Here we will 
consider only minimum phase system [7]. 

Let, 
2.6 1.5 



P^): 
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The interval plants with uncertainty are 
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where, 

^e [2.4, 2.8]; k 2 . 

k 3 <= [1.3, 1.5]; k 4 t 

a x e [58, 66]; a 7 

a 3 e [85, 95]; a 4 



P22 



(l+a^s) 



e [1.4,1.6]; 
; [2.6,2.8] 
e [58,66]; 
e [85,95]; 



The performance specification for tank 1 and tank 2 
are: 1. Tracking specifications: 16 sec. < tn (t) < 27 
sec. and 10% overshoot. 2. Stability specifications: 
GM > 6 dB and PM > 40°. 3. Interaction < 10%. 
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Figure 6: System time domain tracking performance specifications 
for tank- 1. 

The time domain specifications are as indicated in 
figure 6. The time domain specifications are translated 
into the frequency domain. The frequency 
specifications are defined by an upper bound Bu and 
lower bound B L as shown in figure 7. The desired 
control ratios corresponding to upper and lower bound 
are synthesized as the dB difference between upper 
bound (Bu) and lower bound (B L ) increases with 
frequency in high frequency range. The desired spread 
between Bu and B L is achieved 

by augmenting Bu with zero as close as possible to 
origin without significantly affecting the time response. 
The spread can be further increased by augmenting B L 




Figure 7: Frequency domain performance specifications (Bode plot) 
for tank- 1. 

with negative real pole which is as close to origin as 
possible but far enough away not to significantly affect 
the time response. 

The equivalent frequency domain transfer functions 
of upper and lower tracking bounds for tank- 1 are as 
follows: 

0.04^-1-0.04 



H-0.4.S -1-0.04 



B ni 



h-1.528^ 2 -i-0.3757.s- -1-0.02457 



The bounds are generated at different frequencies 
for given stability specification. The disturbance 
specifications are calculated by using equation 15 and 
16. All the bounds are plotted for given stability and 
disturbance specifications. The most dominating 
bounds are selected for controller design at each 
frequency. The bounds on ln and nominal loop 
shaping of lno for tank 1 are as shown in figure 8. 

The corresponding compensator transfer function 
for tank 1 is as follows: 

21.20(5+0.06119) 
gll ~ s(s+ 0.7758)(j+ 8.553) 

The prefilter for tank 1 is as shown in figure 9. The 
corresponding prefilter for tank 1 are as follows: 
0.067 

11 "(5 + 0.067) 



Similar procedure is carried out for tank 2. The 
bounds on I220 and nominal loop shaping of I220 for tank 
2 are as shown in figure 9. The corresponding 
compensator transfer function and the prefilter for tank 
2 are as follows: 

_ 54. 63Q +0.06785) 
g22 ~ sO -1-2.135)0 -1-4. 585) 
The prefilter for tank 2 is 



0.072 



O -l-0. 072) 




Figure 8: The bounds on lno and nominal loop shaping of lno for 
tankl. 




Figure 9: The bounds on I220 and nominal loop shaping of I220 for 
tank 2. 

The SIMULINK block diagram for MIMO system 
is as shown in figure 10. 





Figure 10: The SIMULINK block diagram for MIMO system. 

The closed loop responses of entire plant with step 
input at vi with the designed compensator is within the 
specified upper and lower bounds as shown in figure 
11. Same is true for step input at V2 as shown in figure 
12. 




Figure 1 1 : Closed loop responses yi and y 2 with unit step input at 

Vi. 
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Figurel2: Closed loop responses yi and y 2 with unit step input at v 2 . 

The closed loop response of entire plant with step 
input at Vi and vi with the designed compensator is 
within the specified upper and lower bounds as shown 
in figure 14. 



Figure 14: Closed loop responses yi and y 2 with unit step input at vi 
and v 2 . 

VII CONCLUSION 

In the proposed method, the tracking specifications 
are converted into disturbance specifications, because 
of which the design for MIMO system becomes 
simpler. The EDA method is successfully applied to 
the Quadruple Tank application. The result shows that 
the desired specifications are achieved by both the 
controlled variables. The result also shows that the 
interactions between the loops are within specified 
limits. 



Here, the EDA methods based on QFT methodology 
shows power of QFT. The result proves that QFT 
methods can be used to handle large uncertainties. 
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